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Abstract 

A well-known diffuse interface model for incompressible isothermal mixtures of two 
immiscible fluids consists of the Navier-Stokes system coupled with a convective 
Cahn-Hilliard equation. In some recent contributions the standard Cahn-Hilliard 
equation has been replaced by its nonlocal version. The corresponding system is 
physically more relevant and mathematically more challenging. Indeed, the only 
known results are essentially the existence of a global weak solution and the exis- 
tence of a suitable notion of global attractor for the corresponding dynamical system 
defined without uniqueness. In fact, even in the two-dimensional case, uniqueness 
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of weak solutions is still an open problem. Here we take a step forward in the case 
of regular potentials. First we prove the existence of a (unique) strong solution in 
two dimensions. Then we show that any weak solution regularizes in finite time 
uniformly with respect to bounded sets of initial data. This result allows us to de- 
duce that the global attractor is the union of all the bounded complete trajectories 
which are strong solutions. We also demonstrate that each trajectory converges to 
a single equilibrium, provided that the potential is real analytic and the external 
forces vanish. 

Keywords: Navier-Stokes equations, nonlocal Cahn-Hilliard equations, regular 
potentials, incompressible binary fluids, strong solutions, global attractors, con- 
vergence to equilibrium, Lojasiewicz-Simon inequality. 
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1 Introduction 

The evolution of an incompressible mixture of two immiscible fluids can be described 
through a diffuse interface model (cf., e.g., [HI EQl [221 125] their references). Assuming 
that the temperature variations are negligible, taking the density is equal to one, and 
suppose the viscosity u to be constant, the model H (see [2T]) reduces to the so-called 
Cahn-Hilliard-Navier-Stokes system 

(ft + u ■ Vv? = V ■ (kV/i), 
l^ = -Aip + F'{ip), 

Ut — uAu + (m ■ V)m + Vtt = fj^Vip + h(t), 
div(M) = 0, 

in n X (0, oo), where Q C M*^, d = 2, 3, is a bounded domain. Here u denotes the (average) 
velocity and ip is the difference of the two fluid concentrations. Moreover, k > is the 
mobility coefficient, F is a suitable double well potential density, n the pressure and h a 
given external (non- gradient) force. 

The existing theoretical literature (see, for instance, [H El El [121 (131 [HI [261 129] ) can 
be summarized by saying that all the results known for the Navier-Stokes system can be 
extended to the Cahn-Hilliard-Navier-Stokes one, with some additional technical difficul- 
ties when, for instance, F is a singular (i.e. logarithmic) potential and/or the mobility k 
depends on ip and vanishes at pure phases (cf. [H ISj). However, we recall that the Cahn- 
Hilliard equation has a phenomenological nature (cf. [6]). Instead, a rigorous derivation 
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from a microscopic model yields a nonlocal equation (see [l6l[T7]). In this case the chem- 
ical potential fi has the following form 



where * denotes the convolution product over J : M*^ — > M is a sufficiently smooth 

interaction kernel such that J{x) = J{—x) and a{x) = / k{x — y)dy. Motivated by this 

Jn 

fact, in [7] we have introduced and analyzed the following nonlocal Cahn-Hilliard-Navier- 
Stokes system 

(ft + u ■ V(p = Afi, (1.1) 
f^ = aip-J*ip + F'{ip), (1.2) 
ut - uAu + (n ■ V)n + Vtt = ^Vv? + h{t), (1.3) 
div(M) = 0, (1.4) 

endowed with boundary and initial conditions 

^ = 0, u = ondnx{0,T) (1.5) 
on 

m(0) = uq, (p{0) = (po in Q. (1.6) 

For such a problem we have proven first the existence of a global weak solution satisfying 
an energy inequality (equality in dimension two) for a regular potential F (see [7]). Then 
in PI we have established the existence of a global attractor for the generalized semifiow 
{d = 2) and a trajectory attractor (c? = 3). Similar results have recently been extended 
to singular potentials of logarithmic type (cf. [ID]). However, an important issue has 
been left open: the uniqueness of weak solutions in dimension two. This is well known 
for the standard local models and it suggests that the present model is more difficult 
to handle. The main reason seems to be the poorer regularity of (f which makes the 
capillarity term (i.e. the Korteweg force) difficult to handle (see [7]). Here we are 

not able to address this issue but we come close. More precisely, we prove the existence of 
a (unique) strong solution and the regularization in finite time of any weak solution. The 
latter is uniform with respect to bounded set of initial data so that, as a by-product, we 
deduce that the global attractor we found in [9] is smooth. More precisely, it is the union 
of all the bounded complete trajectories which are strong solutions to fll.ip - fll.6l) . Finally, 
taking advantage of the regularization property, we show that any weak trajectory does 
converge to a unique equilibrium (cf. [121 [221 IM] for nonlocal Cahn-Hilliard equations). 



3 



2 Notation and known results 



We set H := and V := H^iVt). For every / G V^' we denote by / the average of / 

over f2, i.e., / := 1). Here \VL\ is the Lebesgue measure of VL. We assume that dVL 

is smooth enough. 

Then we introduce the Hilbert spaces 

V^:={veV = 0}, ■.= {feV':J= 0}, 

and the operator A : V ^ V , A e C{V, V), defined by 

{Au, v) := / Vn ■ Vv Vn, v G V. 
Jn 

We recall that A maps V onto Vq and the restriction of A to Vb maps Vq onto Vq isomor- 
phically. Further, we denote by A/" : V^' — )■ Vb the inverse map defined by 

AMf = /, V/ e and ^Au = u, Vm G Vo. 

As is well known, for every / G Vq, M f is the unique solution with zero mean value of 
the Neumann problem 

—Am = /, in 
1^ = 0, on dVL. 

In addition, we have 

{AuMf) = {f\u), VuGV, V/G^o', (2.1) 
{fMg) = {g,Mf) = [ V(Ar/) • V{Mg), yf,g G K,'- (2.2) 

We consider the canonical Hilbert spaces for the Navier-Stokes equations with no-slip 
boundary condition (see, e.g., [25] ) 

-L2(n) 



Gd^. := {m G Co~(f^)'^ : div(M) = 0} ' \ Ki„ := G H^iQ)" : div(M) = 0}. 

We denote by || ■ || and (■, ■) the norm and the scalar product on both H and Gdiv, 
respectively. Instead, Vdiv is endowed with the scalar product 



(u, v)va,, = (Vm, Vw), Vm, V eVc 



div ■ 



We shall also need to introduce the Stokes operator S with no-slip boundary condition. 
More precisely, S : D{S) C Gdiv Gdiv is defined as S := —PA with domain -0(5*) = 
H'^{VlY n Vdiv, where P : L'^iVlY Gdiv is the Leray projector. Notice that we have 



{Su, v) = (n, v)v,,^ = (Vn, V^), Vn G D{S), Wv G V 



divy 
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and : Gdiv — ^ Gdiv is a self-adjoint compact operator in Gdiv Thus, according with 
classical results, 5* possesses a sequence of eigenvalues {Xj} with < Ai < A2 < ■ ■ ■ and 
Xj — 00, and a family {wj} C D{S) of eigenfunctions which is orthonormal in Gdiv Let 
us also recall Poincare's inequality 

Ai||m|P < ||Vm||^, Vm g Vdiv 

The trilinear form b which appears in the weak formulation of the Navier-Stokes equa- 
tions is defined as follows 

b{u,v,w) = / {u-\/)v-w, ^u,v,w EVdiv, 
Jq 

and the associated bilinear operator B from Vdiv x Vdiv into V^^^ is defined by 

{B{u, v),w) := b{u, V, w), Vm, v,w e Vdiv 
We shall set B{u, u) := Bu, for all u G Vdiv We recall that we have 

b{u,w,v) = —b{u,v,w), Wu,v,w&Vdiv, (2.3) 
and that the following estimates hold in dimension two 

\b{u,v,w)\ < c||M||^/2||VM||^/2||Vt;||||w||i/2||Vwf/2^ \/u,v,w E Vdiv (2.4) 

\b{u,v,w)\ < c\\u\\^^^\\Vu\\^^^\\Vv\\^/'^\\Sv\\^^^\\wl Vn G Vdiv, v G D{S), w G Gd^v 

(2.5) 

If X is a Banach space and r G M, we shall denote by L^^^ir, 00; X), 1 < p < 00, the 
space of functions / G Lf^^([r, 00); X) that are translation bounded in Lf^^([r, 00); X), i.e. 
such that 

\\f\\lnroo-X)--=^^P / ll/(^)IIV^<00. (2.6) 

We shall use the following lemma. Its simple proof is given below for the reader's 
convenience. 

Lemma 1. Let f G L'^^ij, 00; X) with ft G mij, 00; X), where 1 <pi < 00, 1 < P2 < 00, 
r G M and X is a reflexive Banach space. Then f{t) in X as t —t- 00. 

Proof. We argue by contradiction. Suppose there exist a sequence {tn} with t„ — 00 
and a constant cr > such that ||/(tn)||x > o", for all n. Set r„ := t„ + 1/n. Since 
/ G B^^ij, 00; X) with 1 < pi < 00, then, by possibly extracting a subsequence, for every 
n there exists t'^ G [tri,r„,] such that ||/(t^)||x < cr/2. We therefore get a contradiction, 
since, denoting by p'2 G [1, 00) the conjugate of p2, 

< ^ < ||/(0 - /(t„)|U < r ll/*(^)IUrf^ < ll/*IL-(.,oo;X);;^ ^ 0. 

□ 
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We also report the uniform Gronwall lemma which will be useful in the sequel (see, 
e.g., [27]). 

Lemma 2. Let $ be an absolutely continuous nonnegative function on [t, oo) andui,U2 
two nonnegative locally summable functions on [r, oo) satisfying 

^^(t) < uJi{t)(^{t) + U2{t), for a.e. te [t, oo), (2.7) 

and such that 

rt+l pt+l 

/ Ui{s)ds<ai, i = l,2, / ^{s)ds < a^, (2.8) 
for allt > T, where 01,02,03 are some nonnegative constants. Then 

$(t + 1) < (02 + 03)6"^ , Vt > r. (2.9) 

We now summarize the main results of [7]. They are concerned with the existence 
of dissipative weak solutions and the validity of the energy identity and of a dissipative 
estimate in dimension two. 

The assumptions on J and F are listed below 

(HI) J e lyi'i(M^), J{x) = J(-x), a > a.e. in Q. 
(H2) F G Cf^lO^) and there exists cq > such that 

F"{s) + a{x) > Co, Vs G R, a.e. x e n. 

(H3) F G C^(M) and there exist ci > 0, C2 > and g > such that 

F"{s) + a{x)>ci\s\'^'^ -02, VsgM, a.e. a; G 1^. 

(H4) There exist C3 > 0, C4 > and r G (1,2] such that 

|i^'(s)r<C3|F(s)|+C4, VsGM. 

Remark 1. Assumption J G W^'^{M.'^) can be weakened. Indeed, it can be replaced by 
J G W^'^{Bs), where Bs := {z G M'^ : \z\ < 6} with 6 := diam(f2), or also by (see, e.g., 

m 

sup/ {\J{x-y)\ + \VJ{x-y)\)dy < 00. 
x€n Jn 

The above assumptions allow to prove the following result (see [7]) 
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Theorem 1. Let h E Lf^^{[0, oo); V^^^), uq G Gdiv, (fo E H such that F{(po) E L^{il) and 
suppose that (H1)-(H4) are satisfied. Then, for every given T > 0, there exists a weak 
solution [u,Lp] to fll.ip - fll.6p such that 

u E L°°(0, T; n L2(0, T; ^ E L°°(0, T; L^+^\9)) n L\{), T; V), (2.10) 

WiGL^/3(0,T;V:;j, ^teL'/%0,T;V'), d = 3, (2.11) 

GL2(0,T;\/jJ, d = 2, (2.12) 

E L'^{0,T;V'), d = 2 or d = 3 and q> 1/2, (2.13) 

and satisfying the energy inequality 

S{u{t), ip{t)) + ^ {iy\\Vuf + II V/i||')rfr < Siuo, ^o) + ^ (/^(r), M)rfr, (2.14) 

for every t > 0, where we have set 

£{u{t),ip{t)) = l-\\u{t)f + ] [ [ J{x-y){ip{x,t)-ip{y,t)fdxdy+ I F{^{t)). 
^ ^ Jn Jn Jn 

If d = 2, then any weak solution satisfies the energy identity 

j^£{u,v) + y\\Vuf + II V/if = {h{t),u), (2.15) 

In particular we have u E C([0, oo); G^jt,), E C([0, oo); if) and f^F{(p) E C([0, oo)). 
Furthermore, if d = 2 and h E L^^(0, oo; Vj^^), then any weak solution satisfies also the 
dissipative estimate 

S{u{t),ip{t)) <S{uo,^o)e~''' + F{mo)\n\+K, Vt > 0, (2.16) 

where = (ifQ, 1) and k, K are two positive constants which are independent of the 
initial data, with K depending on Vt, v, J , F and ||/i||r2 cooo-v v 

Remark 2. All the previous results hold for a viscosity v depending on which is 
sufficiently smooth and bounded from above and from below (see [7], cf. also P, ITU]). 
Here we assume v to be constant just to avoid further technicalities in the sequel. 



3 Strong solutions in two dimensions 

In this section we state and prove our main result, namely the existence of a (global) 
strong solution to fll.ip - fll.6p and its uniqueness. More precisely, we have 
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Theorem 2. Let h G Lf ^^{[0, oo); Gdiv) , uq G Vdiv, (po E V H L°°{Q) and suppose that 
(H1)-(H4) are satisfied. Then, for every given T > 0, there exists a weak solution [u,ip] 
such that 

u e L°°(0, T; Vd.,) n L\0, T; H^nf), G L°°(fi x (0, T)) n L°°(0, T; V), (3.1) 



G L2(0, T; G^™), G L2(0, T; if). (3.2) 

Furthermore, suppose in addition that F G C^(M) and that (fo ^ H'^{Q). Then, system 

fll.ip - fll.4p admits a unique strong solution on [0, T] satisfying (13.11) . (13.21) and also 

if e L°°{0,T;W^'P{n)), 2<p<oo, (3.3) 

ipte L^{0,T;H)nL'^{0,T;V). (3.4) 

// J G W^^'-^(M^), we have in addition 

^ e L°°{0,T;H\n)). (3.5) 



Moreover, let v^oi; hi] G V^j^ x H'^{Q) x ^^^^([O, oo); Gdiv), ^ = 1,2, 6e iwo seis o/ data 
and denote by [ui,(pi] the corresponding solutions. Then, there exists a positive constant 
A which is a continuous and increasing function of the norms of the data of two solutions 
and which also depends on T, F, J, Q, v, such that the following continuous dependence 
estimate holds 

\\u2{t)~m{t)\\'' + y2{t)-vi{t)\\l> 

+ [ ||Vn2(r)- Vni(r)fdr+ / || (^2 (r) - <^i(r) frfr 
Jo Jo 

< A(||Mo2 - Moif + 11^02 - ^OlWv^ + \\h2 - /^l |li2(0,T;Gd„)) ' (3-6) 

for every t G [0, T] . 

Remark 3. The regularity properties (I3.1l) - (l3.5p imply that 

u G C([0, oo); Kz„), G C([0, T]-V)n C^([0, T]; H\Q)) 

Actually, we have also G C([0, T]; H^iSl)) for every 5 G [0, 2). Recall that the time con- 
tinuity of the velocity field into Vdiv is a consequence of the fact that u G C^([0, oo); Vdiv) 
and of the following differential identity 

~\\yu\\^ + HSuf + {Bu, Su) = (/iV</p, Su) + (/i, Su), (3.7) 

Zi (JjL 

which is deduced by testing equation (II. 3p by Su. 
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Remark 4. If the condition ipo G L°°(f2) in the first part of Theorem [2] is removed, a 
boundedness estimate for the order parameter ip can still be recovered. In particular, 
it can be proved (see [151 Lemma 2.10]) that for every > there exists a constant 
Cm,to > 0) where m is such that |^q| < m, such that 

sup \\(p{t)\\L^^n) < Cm,to- 
t>2to 

Moreover, fl3.1l) - fl3.4p still hold provided the time interval (0,T) is replaced by (2^0,7"), 
for every T > 2tQ. 

Remark 5. In Theorem |5] condition J G H^^'^(M^) is actually needed to ensure the 
regularity property ip G L°°(0,T; H'^{Q)) only. 

Proof. We shall carry out the proof by providing some formal regularization estimates. 
The argument can be made rigorous by means, e.g., of a Faedo-Galerkin approximation 
technique (see [7] for details). 

We first observe that the property (f G L°°{Q x (0,T)) can be obtained by exploiting 
the same argument used in |31 Theorem 2.1]. Indeed, by multiplying fll.ip by y^ly^l^"^ and 
integrating on Q the resulting equation, the contribution of the convective term vanishes 
due to the in compressibility condition (11. 4p and the proof of [H Theorem 2.1] entails 

sup Mt)\\L^(n)<C, (3.8) 
te{o,T) 

where the constant C depends on the initial conditions, in particular on ||no||, on ||y9o||L°°{n) 
and on T (see [U Estimate (2.28)]). Furthermore, if /i G Lj{,(0, oo; G^ii,) then, thanks to 
the dissipative estimate (I2.16p . we have sup^>o ||</'(^)IU2+2<j(n) < C, the constant C being 
dependent on the initial data and on h only. Hence, due to [H Estimate (2.28)], the 
constant C in (13. 8p does not depend on T. 

As far as the regularity of the velocity u is concerned, notice that, since the Korteweg- 
force term fiVip G L^(0,T; L^(f2)^), then by applying [28| Theorem 3.10], we immediately 
obtain ([XI]) i and (E2])2- 

Henceforth we shall denote by c a positive constant which depends only on J, F and 
Q, while c will denote a positive constant depending on J, F, Q and also on the initial 
conditions uq and ipo (in particular on ||Vmo|| and on ||9?o||l°°(q))- The values of both c 
and c may possibly vary from line to line, even within the same estimate. We shall divide 
the proof into three main steps. 

Step 1. Estimate of (ft in L'^{0,T;H) 

We multiply (II. ip by fj,t in H and get 
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a + F'\v))v'i - (^i, J*^,)+ I (u- Vv)iit + ^^11 V/i||' = 0. (3.9) 



Now, we have 



and 



/ {u-Vip)i2t = / {u-Vip){a(pt- J *^t + F"{(p)ipt) 
Jn Jn 



^ Co II ||2 



c||n||^.||V^f, 



(3.10) 



|(v9f, J *ipt)\ = \{-u ■ Vip + A/i, J * (ft) 
< |(n- V(^, J*v^OI + |(V/i,VJ*(^OI 

<'±y,f + c\\u\\U\Vvf + c\\v^i\\\ 



(3.11) 



Plugging fl3.1UI) and fl3.1ip into fl3.9p . using assumption (H2) and integrating the re- 
sulting estimate in time between and t, we obtain 

^llV/if + I Mdr < i||V^of + j\\\u\\U\^vrdr + c ||V^f rfr, (3.12) 
and on account of the following 



\\Vfir>^\\v^r-cyr, (3.13) 

from (13.121) we are led to the differential inequality 

||V/if < II V/io f + cr +/ rn(r) II V/i(r)frfr, Vte[0,T], (3.14) 

Jo 

where m := c (||m||^2 + l) G ^^^(0, T), for all T > 0. Thus the standard Gronwall lemma 
gives 



V/i G L°°(0,T;i7), VT > 0, 
so that, using also (I3.12p . we infer 

ipeL^iO,T;V), ipteL\0,T;H), VT > 0. 



(3.15) 



(3.16) 



This concludes the proof of (13. ip and (13. 2p . 

Step 2. Estimate of ^pt in L°°{0,T;H) 
We differentiate ( II. ip with respect to time and multiply the resulting identity in H by fit- 
This yields 

/ Vtm+ / l2tUfV^+ I fitU-V^t + W^l^tW" = 0, (3.17) 
Jn Jn Jn 



10 



and, due to fll.4p . we obtain 

/ V^^/i* + II V/itf = / (ftU-Vfit+ / fUfVfit, (3.18) 
Jn Jn Jn 



which entails 



Observe now that 



/ vtm + hvfitf < I {ipy + ifV,). (3.19) 

Jn Jn 



n Jn 
1 d 

2di 
1 d 
2di 



/ ayD^ -{J* ipt, -UfVif-u- Vipt + AyUt) + / F"{ip)iftVtt 
Jn Jn 

a + F"{ip))ipl - (V J * ip^, ut^) - (V J * (ft, uift) 



+ (V J * V/ii) - - / F"\^)^i. (3.20) 
On the other hand we have 

|(VJ*v.„«,^)| < ||VJ|Ui||«,||||v^|Uoc||v.,|| < \\\utr\W,f + c, 
\{y J * ipt,uipt)\ < l|VJ||Li||M||Loo||(pif < ciiMiinHlv^tir) 

|(VJ*v.„V/i,)|<^||V/x,|r + ||VJi|i.||v:.,|r. 
Therefore from fl3.19p we get 

+ + ^l|V/i*ir < c(||n|||. + Hub. + ||n,f + l)||v..ir 

-htf + l I F"\ip)ipl + c. (3.21) 



The integral term containing ip^ can be estimated by means of Gagliardo-Nirenberg in- 
equality in dimension two, that is. 



F''\^)ipi < -c\Mi < -c{y,r + ibtirii vv^tii) < ^w^vtr+ciM' + c. (3.22) 



We now need to estimate V(pt in terms of Vfif In order to do that, let us first control 
V(p in terms of V/i in L^, for every 2 < p < oo. We then take the gradient of /i = 
aip — J * ip + F'{ip), multiply it by Vv?|Vv?P~^ and integrate the resulting identity on Q. 
We get 

[ V<^|V<^|^'-^ ■ V/i = /"(a + F"((p))|V<^|P+ /"((^Va- VJ*<^) ■ V<^|V<^|P-^ 

Jn Jn 

11 



and so, by (H2), we find 



coWVifWl, < ||v</p||i;lv/i|U. + (llVaiUoo + llvJiuoiblUHIVy^lir' 



I LP 



< fllV^r^, + c||V/i||i, + cdlValUo. + wvJh^TMl.- 

We therefore obtain 

||V(/^||lp <c||V/i|Up + c, (3.23) 

witli c depending also on p. We now see tliat tlie L^— norm of V/i can be estimated in 
terms of the L^— norm of ipt- Indeed, using once more the two dimensional Gagliardo- 
Nirenberg inequality, we infer 

< ciiv^f /^i^ii^f < c||v/if /^(ii + M'-'/n 
<c(ii(/.,r-2/p+iiwV(/.ir"2/^+i) 

<c(||^,ir-^/^ + ||n||i;^/-||v^|li;^/^ + i), 

where + q^^ = 1/2 and where we have taken into account fl3.15p and the fact that the 
if ^— norm of /i is equivalent to the L^— norm of — A/i + /i, due to fll.Sp . By fl3.23p we 
therefore deduce the desired estimate 

||V/x|U. <c(l + l|v9i||i-2/P). (3.24) 

We now take the gradient of fit and multiply it in by V(pt- We get 

/ VfXfVipt = I {a + F"{^))\V^t\ 
Jn Jn 



|2 



n 



(Vaift - VJ * ift) ■Vipt+ / F'"{ip)iptVip ■ Vift. (3.25) 
n Jn 

Observe that we have 

F"'{ip)(ptV(p ■ Vipt < c||v?t||L3||Vv2||L6||Vv2t|| 

< c[\\ipt\\ + ytf/''\\v^t\\'^'^ (l + Ibif II V^ill 

<2^\\Vvtr + c\M' + c, (3.26) 
Thus from ()3.25p and (I3.26p and using also (H2), we deduce 

Co 4 4 

12 



so that 

4 



JiV/xtf > \\V^t\\-c\Wt\\-c. (3.27) 





We now go back to (13.211) . By combining (I3.22p and (13.271) we obtain 

\jt I (« + n¥')V? + ll|V/i*ir<«(t)|b*f + c||(/.,r + /3(t) + c, (3.28) 



where a := c(||m||^2 + WuWm + IkilP + 1) and (5 := Hv^Hioo We have a, /3 G /^^(O, T). 

From (I3.28P we can easily infer the desired estimate. Indeed, let us multiply (I3.28P by 
(l + /^(a + F"(<^))(^2)-i 



1 d 



log(l+ I [a + F'\^))^i) < -a{t) + --^ , 



2 +/3(t) + c 

1 



< —a{t) + ^{t) + c\\iptf + c. 

Co 

Integrating this last inequality between and t G (0, T) and using the second of (13.161) 
and the fact that ft{0) G H (since ipo G if^(f2)) we therefore deduce that 

ift G L°°(0, T; i/), VT > 0. (3.29) 

In particular, on account of (I3.23P and (13.240 . we also have 

V/x, Vv9 G L°°(0,T;L^'((])), VT > 0, 2 < p < oo. (3.30) 

Furthermore, by integrating (I3.28P between and t G [0,T] and using (13.270 and (13.290 . 
we also get 

ipteL^{0,T;V). (3.31) 

By comparison in (II. ip we can finally obtain estimates for fi and (f in L°°{0,T; H^{Q)). 
Indeed, we have 

||A/i|| < lly^^ll +c||Vm||||V^||lp, (3.32) 

which implies that A/i G L°°(0, T; ^^(fi)), thanks to (13:29|) and fl330D . Recalling ffT3D 
and the smoothness of dfl, we also have 

fie L'^{0,T;H\n)). (3.33) 
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Apply now the second derivative operator dfj := gf^g^ to (11. 2p . multiply the resulting 
identity by dfjip and integrate on Q. Using the assumption J G W'^'^{M.'^), we get 

/ df^^idlcp= [ (a + F"(^))(9j^)2+ / {d.adjif + djad,ip)dlip 
Jq Jn Jn 

+ / Mja-9jj*¥.)9jy.+ / F"'icp)d.cpdjcpdl 
Jn Jq 

From this identity, by means of (H2) and fl3.30p it is easy to obtain 

\\dif^r>%\dyf--c. (3.34) 

Such estimate together with fl3.33p entail 

ip e L°°{0,T;H\n)). (3.35) 

Step 3. Continuous dependence and uniqueness of strong solutions 

Let us consider two strong solutions zi := [ui,(pi] and Z2 := [u2, ^2] corresponding to initial 
data zqi := [moi,V^oi] and 202 := [mo2,V^02] and to external forces hi and /12, respectively. 
Taking the difference between the variational formulation of fll.ip and f ll.2p written for 
each solution and setting u := U2 — Ui, ip := ip2 — fi, /i := /i2 — /Wi and h := h2 — hi, we 
have 

{ut, v) + i/(Vm, Vv) + h{u2, U2, v) - b{ui,ui, v) = -(v52V/i2, v) + (v?iV/ii, v) + {h, v) 

(3.36) 

{pt, ^) + (V/i, V^/^) = (u2<^2, V^A) - {uipi, V^/^), (3.37) 

for every v G Vdiv and every ip E V. Let us choose v = u and ip = Mp and sum the first 
resulting identity to the second one multiplied by 7, where the positive constant 7 will be 
suitably chosen. After some easy calculations we obtain 

^^Ikf + '^llVnf + b{u2,U2,u) - b{ui,ui,u) + ^^ll'/'llyj 

= -{pVfi2, u) - (<^iV/i, u) + 7(u2, pVAfp) + 7(u, piVAfp) + {h, u). (3.38) 
Notice that 

7(</5, /^) = 7(V5, a(p - J*(p + F\lp2) - F'{pi)) > Co7||v?||^ - li^, J * v) 
> coiyf -l\MvJJ\\v\M > co-f\\pf - \\ipf - c-f^yWy,. (3.39) 



Furthermore, as far as the first two terms on the right hand side of fl3.38p are concerned, 
we have 

|(v'V/i2,«)| < ||¥^||||V/X2|U^lk|U4 < ^l|V«f + c||V/X2||i4||(^f , (3.40) 
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\{ViVii,u)\ = |(/iV</^i,M)| < M\\V^4lA\u\\l^ < -\\Vuf + c\\Vvi\\l,M', (3.41) 
where we have used the bound 

||/i|| = ||a(/. - J * ^ + F'i^,) - < 2||a|U^ M + c||^|| < c||v.||. 

The last two terms on the right hand side of fl3.38p can be estimated as follows 

< yr + ci'MUMi,, (3.42) 

\j{u,^iVAfv)\ < l\\uf + ^y^\\l^yfy,. (3.43) 

Consider the trilinear forms on the left hand side of fl3.38p . By (12.41) we have 

b{u2, U2, u) — b{ui, Ml, m) = b{u, Ml, u) < c\\u\\ II V-uill II Vm|| 
< ^llVuf + c||Vnif ||m|P (3.44) 

Plugging fl339|) - fl3:ii|) into fl3:38|) we get 



' ii«f + iMv,^ + ^11 v^f + ^coyf < c(i + wvifiWl, + II v/i2|ii4)||(/^|p 



1 d / 

+ ci{l\\u2\\l. + Ibillio. +7)||v^|Pv„' + c(7+ llVMif )||Mf + ^M'. (3.45) 
Thanks to fl3.30p . we can now choose 7 = 7^, such that 

r* := Co7* - C(l + ||Vv2l||^oo(o,T;L4(n)) + II^/^2||l°°(0,T;L4(Q))) > 0" 

Hence from (13.451) we deduce 

+ l*Ml^ + ^llVMir + r.ll^f < r^it) (Wuf + 7.11^11^,) + (3.46) 
where 

V := c(||Vmi|P + 7*lk2||?,2 + llv^illioo + 7*) e L\0,T), VT > 0. 
The standard Gronwall lemma then yields 

ll^t)lP + 7*||^(t)ll^j<e^^o''^^^'^(lkoir + 7*ll^o||^j + ^11/^1^^ (3.47) 

for every t G [0, T], where we have set uq := no2 — uqi and ipo := (po2 — (poi- By integrating 
(I3.46P between and t and taking (13.470 into account, we also get 



1 d /, 
2^U'"' 



4 

< 



*||v^f^^^ + 2r. fyfdT 

Jo 

2 II iio ^ 



7*lbofv„' + ^||/^|li.(o,;a,.))(l + 2e^^o.W'^^ / vis)ds), (3.48) 

1 t/ 

for every t G [0,T]. Finally, by combining (I3.47p and (13. 48 p . we obtain (13. 6p . □ 
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Remark 6. It is not difficult to see that the tp component of the strong solution to system 
([Ll])-([L5D satisfies 

^ e C([0,oo);if^(fi)). (3.49) 

Indeed, by combining fl3.17p - fl3.20p and taking into account the regularity properties 
of the strong solution, we can see that [a + F"{Lp))Lpf is absolutely continuous on 
[0, oo). Using (H2) and the fact that tp E C([0, oo); C(fi)) (see Remark[3]) we get ||v?tP G 
C([0,oo)). Now, (ESSD and Ht e Ll^{[0, oo);V) imply that /i G C([0, oo); 1/) and, by 
using f l3.33p again, we also have fi G C^([0, oo); if ^(^2)) so that Afi G C^([0, oo); if). 
Moreover, since u G C([0, oo); L^(f2)) and Vip G Cu,([0, oo); L^(r2)) (cf. Remark [3]), then 
we have u ■ Vip G C^([0, oo); if). Thus from fll.ip we deduce that (ft ^ Cyo{[0, oo); H) 
and, on account of the continuity of t i— )■ ||y9t(t)||, then ipt G C([0, oo); ii). Recall now 
that \/(p G C{[0,oo);H'{nf), for every e G [0,1) (cf. Remark EI). Then, choosing 
e G [1/2, 1), we have G C([0, oo); L^{ny). Thus m ■ G C([0, oo); H) and so ffTT]) 
yields A/i G C([0, oo); H) which entails /i G C([0, oo); ii^(r2)). This and the assumption 
J G W^2,1(k2) ^11^^ deduce ([339]). 

4 Uniform estimates and the global attractor 

In this section we establish some uniform in time regularization estimates by exploiting the 
results proved in the previous section. As a consequence we deduce a regularity property 
for the global attractor of the dynamical system generated by fll.ll) - fll.5p whose existence 
has been shown in [S]. 

Proposition 1. Let h G L^^(0, oo; G^jt,), Uq G Vdi-v, G V H L°°{fl) and suppose that 
(H1)-(H4) are satisfied. Then, the weak solution [u,(p] of Theorem\^ satisfies 

u G L°^(0, oo; Vdi„) n L,yO, oo; H\nf), ip G L°^(l] x (0, oo)) n L°^(0, oo; V), (4.1) 



ut G L^iO, oo; Gdiv), e Lfi,{0, oo; H). (4.2) 

Furthermore, suppose in addition that F G C^(M) and that ipQ G H'^{Q). Then, the unique 
strong solution of Theorem\^ satisfies (14. II) . (14. 2 p and, in addition, 

G i:°°(0,oo;l^^'P(fi)), 2<p<oo, (4.3) 

Gi:°°(0,oo;ii)nL2,(0,oo;V). (4.4) 

If J e iy2,i(M2), we also have 

if e L°°{0,oo;H^{n)). (4.5) 
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Moreover, there exists a constant Ai = Ai(m), depending on m (and on F, J, Q, u), such 
that, for every initial data Zq '■= [uo,ipo] ^ Vdiv x H'^i^), with |(^q| < m, there exists a 
time t* := t*{S{zQ)) > such that the strong solution corresponding to Zq satisfies 

rt+l 

\\Vu{t)\\ + y{t)\\H^n) + \Hs)\\H^n)^ <Ai{m), Vt > f. (4.6) 

Proof. Let us first notice that, setting z{t) := [u{t), (p{t)] and Zq := [uq, (po], by integrating 
the energy identity fl2.15p between t and t + 1 we have 

S{z{t + 1)) + (-llVnf + llV^irjdr < S{z{t)) + Whfdr. (4.7) 

Therefore, using also the dissipative estimate f l2.16p . we get 

/ (2llVM|r + ||V/if jrfr<^(^o)e-''* + F(m)|fi|+K (4.8) 

where the constant K depends on ||/i||L2j^(o,oo;Gdi„) on F, J, Q, u. Notice that the 
initial energy £{zq) can be estimated as 

£{zo) < ^||%f + M||(pof + / i^(<^o), M := sup / \J{x - y)\dy. 

From (14. Sp . setting Ao(m) := F{m)\Q\ + K + 1, we deduce that there exists a time 
to = to{£{zQ)) > 0, given e.g. by to = ^og{S{zo) + c), where S{zq) + c > 1, such that 

[^\\^uf+\\^l^\\'')dr<Ao{m), Vt > tg. (4.9) 

We now establish the uniform in time version of estimates fl3.ip i and fl3.2p i for the 
velocity field. To this aim, notice first that (12. 5p implies (see also |2Sl Lemma 3.8]) 

\\Bu\\ < Cll^f /2||y^||||^^^||l/2^ ^ jj^g^ ^ H"^ {^f H Vdiv- 

Therefore, by splitting the term {Bu, Su) on the left hand side of (13.70 and using the 
estimate above, we get the following differential inequality 

-llVnf + z/||5nf < ^WiiV^W + ^||/if + fxllVnf , (4.10) 

where a{t) := Cj,||m|P|| VmH^. Now, recalling Remark H] (see also the proof of [151 Lemma 
2.10]), the assumption h G L1^^{Q , oo; G div) and the dissipative estimate (I2.16p . we know 
that there exists a constant Co(m) > depending on m, and a time ti = ti{S{zo)) 
depending on S{zo) such that 

sup \\(p(t)\\L^(_n) < Co(m). (4.11) 

t>ti 
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Therefore we have sup^^^ ||yu(t) ||L°°(f7) < Ci(m). Then, using also (14. 9 p and (I3.13p . we get 
||/i(r)Vv9(r)f rfr < C^im), j^'^\{r)dr < Cs{m), (4.12) 

for all t > t2 := max{to,^i}- Therefore, (14.81) and (14.121) allow us to apply Lemma [2] to 
the differential inequality (I4.10p and we deduce that 

\\Vu{t)r < C,{m) := ^(2C2(m) + 211/^11^,^(0 + Ao(m))e^^(-), (4.13) 

for alH > ^3 := ^2 + 1- Furthermore, by integrating (I4.10p between t and t + 1, for t > ^3, 
we obtain 

^ ' ~ ' . ,, . ,,2 



I \Hsm,^^-^ds < C^{m) := (1 + C3(m))C4(m) + -[C2{m) + ||/i||i2^(o,oo;G,„) j ' 

(4.14) 

for all t > ts, where we have also used [281 Lemma 3.7]. Estimates (I4.13P and (I4.14p in 
particular imply (I4.ip i . 

Now, let us write (II. 3p in the form Ut = —Bu — uSu + /xVyj + h and observe that, 
owing to [281 Lemma 3.8] (or (12. 5p ). we have 



t+1 rt+i 



\Bu{s)\\''ds < I \\u{s)f\\Vu{s)\mSu{s)fds < ^(m) := ^C|(m)C5(m), 
for all t > ts, and hence 

utis^ds < CM) ■■= c{cl/\m) + uCim) + C^im) + ||/i||i.^(o,oo;G,,„)) ' (4-15) 

for all t>t3. Note that fHlHD entails (g^Di. 

We are now in a position to get uniform in time regularization estimates for ipt first 
in Iv^ft(r, 00; H) and then in L°°{t, 00; H), for some r > 0. 

Let us note first that, by combining (I3.9l) - (l3.1ip and taking (14.111) into account, we 
obtain the following differential inequality, for all t > ti, 

|||V/i|p + co\M' < {CsMMj,, + c)||V/.|r + C,{m)\\u\\UM'. (4.16) 

Observe that (cf. (gH])) 

t+1 1 

iCsim)\\uis)\\l, + c)ds < Cio(m) := —C,{m)Csim) + c, (4.17) 

Cg{m)\\u{s)\\U\vis)fds < Cn(m) := ^--^C',{m)C5im)Cgim), (4.18) 



18 



for all t > t^. Then, using (14. 9p and fl4.17p . f l4.18p . we can apply the uniform Gronwall 
lemma to fl4.16p in [^3,00) and get 

II V/i(t)f < C^2{m) := {C^^im) + Ao(m))e^^«(-), Vt > U := h + 1. (4.19) 

Now, by integrating fl4.16p between t and t + 1, for t > t4, we also deduce 

rt+l 

Co WvtisWds < Ci3(m) := (1 + Cio{m))Cu{m) + Cii(m), Vt > U. (4.20) 

Estimates 04.191) and fl4.20p imply, in particular, fl4.ip 9 and 04.21) 9. respectively. 
Let us now consider estimate 03.281) . Set 

and notice that, on account of 04. lip , we have 

|||^*(t)f < m < Ci4(m)||^i(t)f , Vt > h. (4.21) 

Then, by arguing as in the previous section and taking 04. lip into account, we easily see 
that 03.28P can be rewritten as follows 

^$(t) + ^11 V^tf < uj{t)^t) + I3{t) + Ci5(m), Vt > ti, (4.22) 

where uj{t) := a{t) + Ci6(m)$(t), and a, /3 the same as in 03.28p . Then, by using 04.2ip . 
OCT]) . Omp and (HIS]), we have 

/ <^{s)ds < Cn(m) := —CrsMCiAm), (4.23) 

Jt Co 

t+l ^ 

u{s)ds < Ci8(m) := c(^-C5(m) + C^im) + Cieim)Cn{m) + 1 j , (4.24) 

t+i 

f3{s)ds < Cig{m) := C^{m)Cr{m), (4.25) 

for all t > By applying once more the uniform Gronwall lemma to f l4.22p in the interval 
[^4, 00), we deduce 

ymi' < C2o{m) := lfc5(m) + Ci7(m) + Ci9(m))e^«(™), (4.26) 

Co V / 

for all t > ts := ^4 + 1. Then, by integrating f l4.22p between t and t + l, for t > t^, and 
using 04.2ip . 04.26P and 03.271) (written with a constant C2i(m) in place of c, for t > ti, 
due to 04. lip ), we also find 

/ \\V^t{s)\\'ds < C22{m) := -(C14C20C18 + Ci9 + C15) + (1 + C',,)C2i, (4.27) 

Jt Co V / 
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for all t > t^, where all Ci depend on m. Observe that estimates (14.261) and (14.271) yield 

gai). 

Furthermore, owing to (I3.23P and (I3.24p . we also have 

||V(/^(t)|Up(f,)2 < C23(m), Vt > ts, 2 < p < oo (4.28) 

Finally, on account of (I3.32p . (I4.26P and (I4.13p . we obtain 

MMh^ < 4 - A/i(t) + fi{t)\\ < C24(m) := c(Ci(m) + Cli\m) + Cl^\m)C23{m)) , 

(4.29) 

for all t > ts, and recalling (I3.34p . provided that J G W^^'^(]R^), we get 

IIv'WIIh^ <C25(m), Vt>t5. (4.30) 
Estimates Km and flCTj) yield (gS]). □ 

Let us now recall the main result about the existence of the global attr actor for weak 
solutions to system (ll.ip -f lT3]) in the autonomous case (cf. Since the weak solutions 
to system (ll.ip - (ll.5p are not known to be unique but the energy identity holds, the 
existence of the global attractor is achieved by using J.M. Ball's approach based on the 
notion of generalized semiflows (cf. [3], to which we refer for the main definitions and 
results). 

We assume that h is time independent, i.e., h G Gdiv, and, for m > fixed, we 
introduce the metric space 

'^m '■= Gdiv X (4-31) 

where 

ym:={veH : F{ip) G L\n), \{ip, 1)| < m}, (4.32) 
The space Xm is endowed with the metric 



d(^2,^l) = ||'«2 - ^^ill + 1^2 - V^ill + / F{(p2)- / F{(p 



1/2 

, VZi, Z2 G Xm, 

Jn Jn 

where Zi := [mi,V5i] and Z2 := [m2,V^2]- 

Suppose that (H1)-(H4) are satisfied and that h G Gdiv Let Qm be the set of all weak 
solutions to system (ll.ip - (ll.6p from [0, oo) to X^ given by Theorem [T] and corresponding 
to all initial data zq G Xm- Then, in [9l Prop. 3 and Thm. 3] it is proved that Qm is a 
generalized semiflow on Xm (i.e., Qm satisfies conditions (H1)-(H4) from [3] in the space 
Xm) which possesses a (unique) global attractor Am- 

Take Zq G Xm and consider a weak solution z := [u, ^p] G C([0, oo); Xm) corresponding 
to Zq. From fl2.14p . written with t = r, we know that for every r > there exists t^. G (0, r] 
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such that z{tr) G Vdiv x V. Thanks to Remark|U we can also assume that ip{tr) G L°°{Q). 
We can therefore write the differential inequality f l4.16p in [tr, oo) and, by integrating 
f l4.16p between t^- and t > t^-, we can see that there exists Sr G (t,-, t] such that ftisr) G H 
and hence (p{sr) G H'^{Q). Summing up, introducing the (complete) metric space 

X^:=V,,,xyl := W e H\n) : 1(^,1)1 <m}, (4.33) 

endowed with the metric 

di{z2,zi) = ||Vm2 - VmiII + \\ip2 - (piWuHn), ^Zi,Z2 G X^, 

then, for every r > 0, there exists G (0, r] such that z{sr) G X^ and starting from 
the time Sr the weak solution corresponding to zq becomes a (unique) strong solution 
z G C {[sr, oo); X^) (cf. Remarks [3] and [6]). Such a solution satisfies the dissipative 
estimate fl4.6p in [st-,oo). Let us consider a bounded in X^ subset B C X^- Choosing 
r = 1 for every Zq & B, then every weak solution z starting from Zq & B becomes (at a 
certain time si G (0, 1] depending on zq and on the weak solution considered from zq) a 
strong solution satisfying (14.61) in [l,oo). We therefore deduce that there exists a time 
t* = t*{B) > 1 such that 

z{t) e Bi{Ai{m)), yt>t*, (4.34) 

where Si(Ai(m)) is the closed ball in X^ given by 

B^{A^{m)) := {w e X^ : d^{w,0) < Ai(m)}. 

This fact immediately implies that Am C Bi. Indeed, we have dist;t'i^(T(t)^m, i3i) = 

^1 

distx^{Am, Bi) = 0, which implies Am C Bi = Bi. We recall that the multivalued 
evolution map T{t) is defined, for every t > and every subset E C Xm, as (cf. [3]) 

T(t)E := {z{t) : z e Gm, z{0) G E}. (4.35) 

Summing up we have just proven the following regularity result for the global attractor 

Theorem 3. Let (H1)-(H4) be satisfied and assume that h G Gdiv is independent of 
time. Then the global attractor Am of the generalized semiflow Qm associated with system 
ffrT]) - ffT3]) is such that 

AmCBi{Ki{m)). 

Thus the global attractor is the union of all the bounded complete trajectories which 
are strong solutions to fll.ip - fll.6l) . 
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5 Convergence to equilibria 



In this section we shall prove that every weak solution to system fll.ip - fll.6p converges to 
a stationary solution as t — ?■ oo, provided that F is real analytic and h = 0. 

Let us first consider the set of all stationary solutions z^o to system f ll.ip - fll.Sp . namely 
the set of pairs z^o := [0,v?oo] G '^m (for some m > 0), where (poo solves the integral 
equation 

a(poo- J *^oo + F'{ipoo) = jJ^OO, (5.1) 
with some constant /ioo G M given necessarily by /ioo = F'^ipoo)- Therefore we introduce 

Sm = ^Zoo = [O,ipoo] ■ (poo&H, F{(poo) e L^{Q), < m, 

av^oo - * V^oo + -^'('/'oo) - -^'(V'oo) = a.e. in (]|. (5.2) 

We point out that, by using an easy iteration argument from (15. ip . on account that F' has 
polynomial growth, we can deduce that ipoo G L°°(f2). The structure of the stationary set 
is rather complicated. In particular, such a set may be a continuum (see [8] for an example 
and [T^ where the author proves the existence of solutions ipoo to fl5.7p with ip^ = in 
cylindrical bounded domains). It is also worth observing that to every stationary solution 
2^00 = [0, v^oo] there corresponds a stationary pressure tToo given by tToo = F'{(poo)(poD + c, 
where c G M is an arbitrary constant (cf. (11. 3p ). 

We begin with the following preliminary but crucial result. 

Lemma 3. Assume that (H1)-(H4) are satisfied. Take zq G and let z G C([0, oo); Xm) 
be a weak solution corresponding to Zq. Then, we have 

^ uj{z) C (5.3) 

and 

u{t) —7-0 in Gdiv, as t — )■ oo. (5.4) 

Furthermore, there exists a time t* = t*{zo) depending on Zo such that the trajectory 
[Jt>t*{^{'t)} precompact in Xm- 

Proof. From (I2.14p . by letting t oo, we obtain that 

M G L2(0,oo;1/rfi„). (5.5) 
On the other hand, from (II. Sp . written as Ut = —Bu — vSu + fJ^Vip, we get 

< i^||Vu|| + c||u||||Vm|| + ||(/?||lcx=(q)||V/x||. 
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Now, f l2.14p also implies that u G L°°(0, oo; Gdiv) and that VyU G iv^(0, oo; H). Hence, on 
account of (14. lip as well, from the previous estimate we infer that 

UteL\r,^-XJ. (5.6) 

for some r > 0. From f l5.5p and f l5.6p we deduce ( 15 ■4p . Let us now take z G uj{zq) arbitrary, 
with z := [n, ip\. Then, there exists a sequence with t„ — oo such that u{tn) — > u in 
and (^{tn) if in H. We get u = and, up to a subsequence, 

l^{tn) -> yW, a.e. in (5.7) 

where yU := a(f — J*(f + F'{ip). By integrating (I4.22p between t and t + 1 we easily deduce 
that V/it G Lj^(r, oo; iJ) for some r > 0. Since we also have V/i G L^(0,oo;iJ), then 
Lemma [1] yields 

V/i(t) ^0 in H, as t ^ oo. (5.8) 

From (15. 7p and (15. 8p we easily deduce that /i=const almost everywhere in Q, where 
the constant is necessarily given by F'ijp). Therefore z = [m, = [0,"^ G Em (note 
that F{ip) G L^{VL) is ensured by Fatou's lemma), and (15. 3p is proven. Finally, the 
precompactness of the trajectory is an immediate consequence of (I4.34p . □ 

Remark 7. Lemma |3] yields in particular an existence result for equation (15. ip . 

We now recall the generalized Lojasiewicz- Simon inequality established in [TT] which 
is the main tool for proving our convergence result. 

Let V and W be Banach spaces embedded into a Hilbert space H and its dual if', 
respectively, with dense and continuous injections. Assume that the restriction of the 
Riesz map R G C{H, H') to V is an isomorphism from V onto W = R{V). Moreover, let 
H = Hq + Hi, where Hi G V is a. finite-dimensional subspace and Hq is its orthogonal 
complement in H. Introduce the subspace of H' 

H^ ■={geH' : {g,ip)=0 for all ipeHo}. 

Then let 

where the functionals Qi and Q2 satisfy the following conditions 

• Qi '■ U C y — 7- M is Frechet differentiable on an open set U such that the Frechet 
derivative DQi : U W is a. real analytic operator which satisfies 

(D^l(V92) - Dgi{Lpi),Lp2 - ifl) > ai\\Lp2 - ^i\\h, (5.9) 



23 



DGi{V2) - DGi{<Pi)\\h' < a2\\<^2 - <^i\\h, 



(5.10) 



for all V9i,(/?2 € U and for some constants ai,a2 > 0. Furthermore, the second 
Frechet derivative D'^Qi{ip) G C{V,W) is assumed to be an isomorphism for all 

• Q2 : H ^ assumed to be in the form 

where /C G C{H, H') is a self-adjoint compact operator such that its restriction to 
y is a compact operator in £(V, W) and / G W , p G M are given. 

The inequality we need is given by 

Lemma 4 (p3]). Let the previous assumptions he satisfied for the spaces V, W, H, H' and 
for the functional T . Let [v^oc/^oo] G ?7 x satisfy DJ^{ip^) = fi^. Then, there exist 
0", A > and 9 G (0, 1/2] such that the following inequality holds 

I J-(^) - J-((poo)r"' < A inf { \\DF{^) - pIIh', /i G i^o}, (5.11) 

for all if eU satisfying ip — ipoo £ Hq and \\ip — ifooWu < o"- 

We can now state the main result of this section. 

Theorem 4. Assume that (H1)-(H4) are satisfied with F real analytic. Take zq G Xm 
and let z G C([0, 00); Af^) he a weak solution corresponding to zq. Then, there exists 
Zoo '■= [0, V'oo] G £m with Tp^ = TpQ such that 

z{t) Zoo in Xm, as t ^ 00. (5-12) 

Moreover, there exist some constants 7 > 0, 9 G [0, 1/2) and a time t > which depend 
on Zq and Zoo ( and on the weak solution z originated from zq ) such that 

\Ht)\\v:,^ + Mt) - ^ooWv < it-T^, vt > t. (5.13) 

Proof. Our aim is to prove that ipt G L^{t, 00; V'), for some r > 0. This, together with 
f l5.4p and with the precompactness of the trajectory in Gdiv x H, will allow to deduce 
the convergence in Gdiv x i7 of a whole trajectory z = [u, p] originating from an initial 
datum Zq = [uq, ipo] G Xm. to a stationary solution z^o G Sm with Ip^ = TpQ. Observe that 
if z : [0, 00) — 7- Xm is a weak solution, then the convergence condition z{t) — z^o in Xm is 
equivalent to the condition z(t) — )■ z^o in Gdiv x H, since the convergence J^F{ip{t)) — )■ 
F{ipoo) is ensured by (14. lip and Lebesgue's dominated convergence theorem. 
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The key point is the apphcation of Lemma H] to a suitable functional T which is, in our 
case, the energy functional E associated with the component of the solution, namely, 

E{^) = \\^a^f-\{ip,J^^)-r I F{^). (5.14) 

More precisely, we set (cf. Lemma Hj) 

H ■= H' = L'^in), Ho:= {tp e H -.^ = 0}, iJ° = {^/^ = const}, 
V = L^{n), W := R{V), WfWw ■■= \\R-'f\\v, 

:= ^ + ^a^'), U = U^:={^eV: \i;{x)\ < Co(m), a.e. x G Q}, 

]C{i)) := - J*ip, / = p = 0, (5.15) 

where the positive constant Co(m) is the same as in (14. lip . 

All the assumptions of Lemma H] are fulfilled. Indeed, Qi is Frechet differentiable on 
the whole V with DQi{(f) G W, for all (f & V given by 

{Dg^{^),h)= [ {F'{ip) + aif)h, \/heV. 
Jn 

Furthermore, DQi is a real analytic operator, since F is assumed real analytic, and we 
have 

{DQi{(p2) - Dgi{ipi),(p2 - Vi) = / {F"{ri(p2 + {I - ri)ipi) + a)\ip2 - ^il"^ 

Jn 

thanks to (H2), where r] = rj{x) G (0, 1). Hence (15.91) is satisfied (with ai = cq). As far as 
(I5.10p is concerned, observe that DQi is locally Lipschitz from V to H'. Indeed, we have 



\\Dg,iip2) - Dg^iip^)\\H' < \\F'iip2) - F'(^i)|| + aj\ip2 - ipi\\ < T^\\ip2 - ^if , 

for all fi,(p2 G Um, which yields (I5.10p (with a2 = F^). Moreover, the second Frechet 
derivative is given by 

{D^g,{if)h,,h2)= I {F"{^) + a)h,h2, ^h,h2eV, 
Jn 

for all (yj G V. Hence D'^Qi{ip) G C{V,W) is an isomorphism for all cp G Um- Finally, 
thanks to (HI), the convolution operator /C is compact from H to H and also from V 
to W (due to the compact embedding W^'°°{Q) "^M- C{Q)). The Frechet derivative of 
= E is given by 

DE{<p) = F'{<p) + a<p - J *^ = fi, (5.16) 
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and we have that [v^ocyUoo] e Um x satisfy DE{(poo) = /ioo iff 2:00 := [0,v?oo] G Sm 
with (poo G t/m and /ioo = -^'(v^oo)- Therefore, taking [(/?oo,Aioo] G f/m x Hq such that 
DE{(pac) = /ioo 5 Lemma m entails the existence of a, A > and 9 G (0, 1/2] such that 

\E{ip) -^(v^oo)^"^ < Ainf {||/t-/i||, /i = const} = A||//-/l|| < Acp||V/i||, (5.17) 

for all ip G t/m satisfying ^ = ip^ (i.e. — (poo G i/o) and ||v? — V'ooIIh < o") where Cp is 
the Poincare-Wirtinger constant. 

Now, let zq G Xm. and z be a weak solution corresponding to zq. Take Zoo G i^(^) and 
let {tn} be a sequence such that t„ — j- 00 and ^(^„) — )■ 2:00 in A'^. Consider the function 

$(t) :=^(^(t))-^(^oo). 

We have 

$'(t) = -z/||Vn|p - ||V/i|p < -c^(||Vn|| + < 0, for a.a. t > 0, (5.18) 

where Cj^ = min{l, i^}/2. Since $(tn) — ^ and $ is non-increasing in (0, 00), then $(t) — > 
0, as t ^ 00 and $ > 0. Now, due to ([53D and to fl5T7D (notice that 2(1 - 9) > 1), we 
have 

^'-'(t) = (i||n(t) f + E{^{t)) - E(^oo)) 

< \\uit)r^'~'^ + \Ei^it))-Ei^^)\'-' 

< ca(||Vm|| + ||V/i||), (5.19) 

for all t > to, for some to > 0, provided that ||</?(t)— <y5oo|| < o", where c\ = max{l/-\/Ai, Acp}. 
Therefore, by combining f l5.18p and (15.1 9p we get 

- j^^%t) = -9^'~\tW{t) > ^-^{\\Vum + ||V/i(t)||), (5.20) 

provided that (p(t) G Um with \\(p{t) — y9oo|| < cr and Tp(t) = Tp^ = TpQ. By means of 
a classical argument, together with equations (11.11) and (11.21) . we can now deduce that 
ipt G L^{t, oo]V'). Indeed, for every 6 G (0,1) there exists N = Ns such that for all 
n > Ns we have [[^(tn)!! < 6 and ||v5(tn) — fooW < 5- Set 

r = r(5):=sup{t>tjv: lk(s)|| <1, yis) - ip^\\ < a, VsG[tjv,t]}. (5.21) 

Then, estimate (I5.20p holds for all t G [tN,t*]. By integrating it between In and t* and 
possibly choosing a larger N we have 

J^' (llVii(r)ll + ||V/i(r)||)rfr < ^^^'M < S. (5.22) 
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We now claim that there exists 5* > such that t*(5*) = oo. Indeed, suppose this is not 
true, i.e. t*{6) < oo for all 6 > 0. Then, we have 



\\utir)\\viJr < / (^HlVn(r)ll +c||n(r)||||Vn(r)|| + ||^(r) || V^(r) || jdr 

< h f (\\Wu{t)\\ + \\W^l{r)\\\dr < h6, (5.23) 



where bi = max {i^ + cAi(m) / Co(m)}, and where Ns is assumed large enough, i.e., 
such that ^AT^ > ti{S{zQ)) (see (14. lip ). Furthermore, we have 

* yt{r)\\v'dr< f (||VMr)|| + ||^(r)|U^||n(r)||)dr 

ft* 

<b2 (llVM(r)ll + ||V/x(r)||)rfr < b^S, (5.24) 
where 62 = max |l, Co(m)/A/Ai}. Therefore, we deduce 

MnWv;^^ < HtN)\\vi^^ + r \\utir)\\viJr < b^S, (5.25) 

JtN 

\\^{t*) - ifiocWv < \\^{tN) - ^ocWv + \\^t{T)\\v'dT <b46, (5.26) 

JtN 

where 63 = 1/y/Xi + bi and 64 = 1 + 62- Let us now take a sequence such that 5„ — )■ 0. 
Then, from definition fl5.2ip . for every n at least one of the following two conditions holds 

||n(r((5„))|| = 1, Mt*{Sn)) - VooW = a. (5.27) 

By possibly extracting a subsequence we have, e.g., \\ip{t*{6n)) — V^oo|| = o'. Writing (15.261) 
with 6 = 6n and taking into account the precompactness of the trajectory in Gdiv x H we 
get a contradiction. Thus, for some 5* > we have (setting t := t^gj 

j_ (||VM(r)|| + ||V/i(T)||)c^r<5, <oo, (5.28) 

so that 

u ^ L^(t,oo]Vdiv): Vfi e L\t, 00; H). (5.29) 
This implies that (ft £ L^(i, 00; V), due (I4.ip 9 and to the estimate 

W^tWv < l|V/i|| + C||v2||y|| Vm||. 
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By using the precompactness of the trajectory in Gdiv x H again, we deduce that if{t) — )■ 
(yjoo in i7 as t — )■ oo. Therefore we have z{t) — )■ z^o in Xm as t — )■ oo. We now provide an 
estimate for the convergence rate in V^'^^ x V . Indeed, from flS.lSp and fl5.19p we deduce 

which yields by integration 

< $(0){l + 65$^-2^(0)t}"^, Vt>t, (5.30) 

where 65 = Cy{l — 26')/c|. On the other hand, by integrating fl5.20l) from t > t to 00 we 

get 

j'^ (||V«(r)|| + l|V/i(r)||)rfr = Vt > t. (5.31) 

Finally, we obtain 

\Ht)\\v^^^ < r \\ut{T)\\y,jT < h r (llVn(r)ll + ||V/i(r)||)c/r, (5.32) 



t Jt 

00 POO 



Mt)-ip^\\v'< \Wt{r)\\v'dT<h (||Vn(r)|| + ||VMr)||)cir. (5.33) 
By combining fl5.30p - fl5.33l) we deduce the convergence rate estimate fl5.13p with 7 = 

Remark 8. By using standard interpolation inequalities one can deduce from f l5.13p 
convergence rate estimates in stronger norms. Of course, the convergence exponent dete- 
riorates. 
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